Incomplete electronic relaxation and population upconversion in quantum dots 
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We calculate the electronic relaxation of a single electron in a quantum dot with two electronic 
orbital states, and with the electronic coupling to the longitudinal optical modes of the lattice 
vibrations included in the multiple scattering approximation to the electronic self-energy. This 
model shows that there is an incomplete electronic relaxation from the excited state to the ground 
state and also the incomplete increase of the population of the excited state after resonantly exciting 
the ground state. These theoretical results are compared with very similar recent findings reported 
in experimental papers. 
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I. INTRODUCTION 

Relaxation of the electronic energy in individual quan- 
tum dots and the question of the phonon bottleneck still 
remain to be of considerable current interesliiSi^. The 
Frohlich's mechanism of the electron-LO-phonon cou- 
pling in the polar semiconductors is regarded as very 
important^.. The electronic self-energy due to this elec- 
tronic interaction was studied by Inoshita and Sakaki^, 
who found that the spectral features of the electronic 
spectral density have very narrow width. The depen- 
dence of these very narrow spectral peaks on the energy 
parameter E was shown to have the shape of the E~^l'^ , 
in the limit of zero temperature^. 

The relaxation of the electron energy in quantum dots 
was studied lateB£*2*24i& in the approximation of multi- 
ple electronic scattering on optical phonons expressed in 
the self-consistent Born approximation. It was theoreti- 
cally shown that despite the fact that the spectral peaks 
of the electronic spectral density have zero width in the 
sense of full width at half maximum, the relaxation rate 
of an electron from the excited state in the quantum dot 
remains nonzero in the limit of zero temperature of the 
latticeZ*2i. The relaxation rate obtained in this way was 
shown to have the order of magnitude of l/ps2iS, in agree- 
ment with broad experimental data. 

In the same approximation to the electronic self-energy 
it was demonstrated that the model gives the electronic 
relaxation even in cases when the electronic inter-level en- 
ergy separation is not equal to the optical phonon enery. 
The dependence of the relaxation rate on the detuning, 
between the electronic energy level separation and the en- 
ergy of the optical phonon, has a resonance structureii. 
The existence of this resonance structure is in accord with 
the well-known experimental rule for detecting the pres- 
ence of quantum dots in quantum dot samples with help 
of optical method&4^ basing on the detection of the reso- 
nance structure of the photoluminescence excitation sig- 
nal with peaks separated by optical phonon energy Ej^o- 

The important role of Frohlich's coupling has also been 
supported by the recent interpretation of the optical line 
broadening and of the continuous background observed in 
photoluminescence excitation spectra measured on self- 



assembled quantum dot samples at high levels of the elec- 
tronic excitation^^. 

Treating the electronic relaxation in the self-consistent 
Born approximation to the self-energy means that the 
electron makes the irreversible transition between two 
electronic states which are renormalized with respect to 
the electron-phonon coupling, with the electronic self- 
energy taken in the self-consistent Born approximation. 
This approximation takes into account multiple scatter- 
ing processes of electron on the system of optical phonons 
and introduces in this way the multiphonon states into 
the relaxation rate formula and into the formula for the 
self-energy. The coherent multiphonon states created 
in this way are analogous to the classical macroscopic 
oscillations^^ of the lattice and the irreversible process of 
the electronic relaxation in this system can be viewed as 
a quantum transition in a potential dependent on timei^. 
A nonconservation of energy in such a system then may 
be manifested as an absence of the phonon bottleneck in 
quantum dots. Because of the rapidity of the electronic 
relaxation in quantum dots, this property of the elec- 
tronic relaxation between the electronic states in the dot 
restricts considerably the possible use of quantum dots 
in constructing quantum bits using the electronic orbital 
states in quantum dots as qubit statesi^ii^. 

To our knowledge, in the theory the electronic relax- 
ation has so far been studied with an emphasis to the 
relaxation rate from the excited state only, calculated 
at the innitial stage of the relaxation, which was under- 
stood to be such instant of time, at which the electron 
is created at this stateSiS. The time dependence of the 
electronic relaxation in quantum dots has not yet been 
given much attention in theory. In this paper we present 
the numerical evaluation of the time dependence of the 
relaxation. We point out several features which are seen 
from the numerical data, concerning the final stage of 
the relaxation. In particular, we demonstrate theoreti- 
cally the effect of the rapid decrease of the population of 
the ground state energy level in the quantum dot and an 
appearence of electronic population in the excited state, 
after resonantly exciting the electron to the ground state 
(an upconversion). The effect of the incomplete depop- 
ulation of the electronic excited state in quantum dot. 
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at the end of the electronic relaxation from the higher- 
energy state and at low temperatures, is also shown. We 
point out possible connection of this theoretical effect to 
the recent experimentsiii2*iS. 

II. THE TRANSPORT EQUATION 

Similarly, as it was done in the previous works, we re- 
strict ourselves to a simple model of the electronic and 
phonon structure of the quantum dot. First of all, we 
ignore completely the presence of the holes in the va- 
lence band states of the dot. Although the holes may be 
simply assumed as being in thermal equilibrium with the 
surroundingSS, we simplify our model assuming the holes 
to be heavy enough, so that they do not influence consid- 
erably the motion of the system under consideration and 
contribute only to the static potential of the quantum 
dot. 

We consider the presence of a single electron only. 
The quantum dot electronic orbital states, and the corre- 
sponding energies, are approximated by the eigenstates 
of the infinitely deep cubic quantum dot, with the lateral 
size d and with the material properties of GaAs crystal. 
We consider again the two-energy level modeliS, assum- 
ing that the electronic ground state ipQ (r) , with the index 
n = and with the energy Eq, is given by the electronic 
ground state in the infinitly deep cubic quantum dot. 
The excited state, with V'l (r) being the wave function of 
the state n = 1 with the energy Ei , is chosen to be one of 
the triple degenerate lowest energy excited states in this 
cubic dot. We put i?o = in this work. 

The electron is assumed to interact with the bulk 
modes of the longitudinal optical phonons of lattice vi- 
brations, neglecting in this way any influence which 
the interfaces in the quantum dot may have on the 
structure of the optical phonons. It was shown ear- 
lier that the details of the optical phonons can be ne- 
glected in heterostructures with not too small character- 
istic dimensions^. We confine ourselves to this simple 
model, neglecting all other interaction, like the electronic 
coupling to the acoustic phonons, and other electronic in- 



teractions. Neglecting the electronic spin, the Hamilto- 
nian of the electron-phonon system under consideration 
readsifii^S 

H = Ho + H,, (1) 

where Hq is Hamiltonian of free electrons and free optical 
phonons, 

n=0,l q 

Cn is annihilation operator of electron in state n, 6q is 
annihilation operator of optical phonon with wave vector 
q. The electron-phonon interaction operator (Frohlich's 
coupling) is 

m,n— 0,1 q 

where the summation over q extends over the phonon 
wavevectors in the first Brillouine zone. In the latter op- 
erator the coupling constant is Aq — {—ie/q)[Ei^o{'^^ ~ 
Kg ^)]-^/^(2eo^)^^^^, where and kq are, respectively, 
high-frequency and static dielectric constants, eo is per- 
mittivity of free space, — e is the electronic charge, 
g =1 q I, and V is volume of the sample. $ is the form- 
factor, 

*(n,m,q)= j d'vi^l{v)e''''iJm{v), (4) 

which modifies the original bulk form^i of the electron- 
phonon operator to the case of quantum dot. 

The relaxation rate dNi/dt can be obtained using 
the nonequilibrium Green's functions method—. As- 
suming the self-consistent Born approximation for the 
electronic self-energy, the instantaneous collisions ap- 
proximation, the diagonal approximation for the single- 
particle Green's functions and distribution functions, and 
the Kadanoff-Baym ansatz-- with the use of the equi- 
librium spectral densities for the interacting system, we 
finally get: 



dNi 
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(5) 
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We do not give here the details of the derivation of 
the relaxation rate formula, referring the reader to 
referencesSiiii In formula ^ crn{E) is spectral density 
of the state n calculated in the self-consistent Born ap- 
proximation, E being the energy variable, v^o is Bose- 



Einstein distribution, giving the population of the optical 
phonon modes at the given temperature T of the lattice. 
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The constant is given by the definition 
q 

The spectral densities, obeying the condition 
J^^an{E)dE — 1, are determined from the corre- 
sponding Green's functions. The retarded electronic 
self-energy M„(_E) of the state n is obtained, in the self- 
consistent Born approximation, from the self-consistent 
equation 

Mn{E) = ^a„„ 

m 

^ / l-N^ + vlo 

"^yE-E,,,- Elo - M„,iE - Elo) + «0+ 

E-E^ + Elo - {E + Elo)+^Q+ ) ' 

In the formula for the relaxation rate the spec- 
tral densities depend on the electronic density matrix by 
means of the equation ^ for the self-energy. This de- 
pendence of the relaxation rate on the state of the elec- 
tronic system is not easily seen in a quantitative manner 
from the formula (O and will be therefore demonstrated 
numerically below. Let us remark, that in the case of 
approximating the electronic spectral densities by delta- 
functions, which would be a suitable approximation in 
the case of very weak electron-phonon interaction, the in- 
tegrals on the right hand side of (|SJl would be performed 
and the dependence of the relaxation rate on the popu- 
lation of the unperturbed electronic states would be be 
simple. This approximation would give nonzero relax- 
ation rate only when the inter-level energy separation 
equals the optical phonon energy. We shall therefore not 
pay any further attention to this simple approximation. 

III. ELECTRONIC RELAXATION 

In the present approximation, in which we confine our- 
selves to the diagonal terms of the electronic density 
matrix, given by Nn, n — 0,1, the electronic density 
matrix is determined by one of the terms -/V„, because 
Nq -\- Ni — 1. This fact allows us to plot the electronic 
relaxation rate —dNi/dt simply as a function of Nq. The 
dependence of the relaxation rate on A'o is shown in 
the Fig. n The plot is shown for three values of the 
quantum dot size, in which the corresponding inter-level 
separation, or the detuning, is set up to be such, that 
El = ^Elo, with ^ being 0.7, 1 and 1.7. The value of 
^ = 1 corresponds to the resonance between Ei and E^o- 
The relaxation rate is not constant as a function of A^o- 
Also, the relaxation rate depends on the magnitude of 
the detuning ^. The relaxation rate changes sign at some 
values of iVg, which are different from and 1. Another 
numerical result says that the dependence of the rate on 
the population is monotoneous. We see in Fig. ^that at 
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FIG. 1: Electronic relaxation rate —dNi/dt plotted for three 
values of the detuning, ^ = 1, 0.7 and 1.7, as a function of 
the population A'o of the lower-energy state n — 0. Lattice 
temperature is T=10 K. 



all three cases of the detuning ^ the relaxation rate be- 
comes zero at the population, which is about between 0.1 
and 0.2 at the lattice temperature of 10 K. Realizing that 
the optical phonon energy in GaAs is about 36.2 meV, 
the electronic inter-level energy separations considered 
in this figure are of the same order of magnitude. The 
value of Nq in Fig. ^ at which the rate is zero and the 
relaxation stops, is therefore not expected to be simply 
due to the equilibrium effect of the lattice temperature 
electronic level population. 

Fig.^shows that the final state of the electronic relax- 
ation from the state with the population A'^i = 1, at low 
temperatures, cannot be the state with iVi = 0, but it 
should be another state of the sytem, in which the pop- 
ulation of the ground state n = is different from one. 
In other words, the electron does not relax to the ground 
state even in the limit of low temperatures. In Fig. |2|we 
plot the final population iVi / of the excited state, which 
is defined as such a value of the population A'^i , at which 
the relaxation rate becomes zero. The quantity Nif is 
plotted as a function of the detuning and lattice temper- 
ature. We expect that due to the population dependence 
of the relaxation rate, shown in Fig. ^ the population of 
the state n = 1 will be always either decreasing at i > 0, 
if we start with A^i = 1 at t = 0, or increasing, if we 
start with the population A^i = at < = 0. The popula- 
tion of the states will be developing towards the limiting 
population Nif plotted in Fig.|21 

The final value A^i / of iVi depends both on lattice tem- 
perature and on the detuning ^. The quantity Nif ap- 
pears to increase with increasing the lattice temperature, 
which is expected. The dependence of Nif on the detun- 
ing ^ displays a certain modulation. At the integer values 
of ^, which correspond to the resonance between the ex- 
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FIG. 2: The final population Nif, at which the relaxation rate 
is zero, as a function of detuning ^ and lattice temperature 
T. 
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FIG. 4: Time evolution of the population A'^i, of the higher- 
energy state n = 1, at the detuning ^=0.7. 
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FIG. 3: Time dependence of population iVi of the higher- piG. 5: Time evolution of the population A^i, of the higher- 
energy electronic state n = 1 for the detuning ^=1, (Ei = energy state n = 1, at the detuning ^=1.7. 
Elo), at lattice temperature T—10 K. 



citation energy Ei and a multiple of the LO phonon en- 
ergy, the value of A^i / shows a minimum, which is rather 
pronounced especially in the limit of low temperature 
and near = 1- At higher lattice temperatures the lat- 
ter resonance structure appears to weaken. The present 
two-level model shows that the limiting value iVi / is es- 
pecially large, achieving the magnitude of about 0.4 to 
0.5, when the detuning is either much larger, or much 
lower than 1. 

Studying the time dependence of the electronic distri- 
bution in the process of the electronic relaxation, let us 
first assume that the electron is prepared in the state 
71 = 1 at t = 0. The numerical result is displayed in Fig. 
13 On the time scale of picoseconds, the population of 
the excited state decreases with time, reaching the limit- 
ing value of the electronic population, in agreement with 



the data shown in Fig. 12 

It is very interesting to observe the time development 
of the electronic distribution in the present two-level 
model in the particular case when the electron is pre- 
pared in the ground state n — 0, namely when the pop- 
ulation Ni = is prepared at time t = 0. This innitial 
condition may correspond to the experiments, in which 
time resolved depopulation of the energy levels is mea- 
sured, including the time dependence of the optical prop- 
erties of quantum dots after exciting the quantum dots 
to the lowest energy state of the electronic subsystem 
in a single quantum dot^iiS. The time dependence of 
the population A^i corresponding to this initial condition 
is displayed in Figs. 01 and |5l in which this time depen- 
dence is shown for two values of the detuning ^. From 
these figures it is seen that on the time scale of picosec- 
onds the population of the electronic excited state can 
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increases to about 0.2 to 0.3, which means that the pop- 
ulation of the lower energy state can decrease to 0.8 to 
0.7. This happens in the limit of low temperatures, at 
which the thermal population of the higher-energy states 
would be unexpected. From the difference between Figs. 
0]and|Slwe see that the characteristic time of this pop- 
ulation upconversion slightly varies with detuning and 
temperature. In agreement with the data shown in the 
Figs, m the results shown in Figs. 01 and [S] demonstrate 
that at temperatures below about 100 K the rapidity of 
the population changes does not depend significantly on 
temperature. 

IV. DISCUSSION 

The present theoretical results can find a counter- 
part in the recent experiments in which the time re- 
solved response of quantum dots to short laser pulse is 
studiediii^. The experimental experience shows that 
when a quantum dot is excited with a short laser pulse 
to the lowest-energy state, the subsequent relaxation pro- 
cess may lead to the appearence of the electronic popu- 
lation of electronic states with higher energjii^. Another 
effect found in experiment is: the system of quantum 
dot, when prepared in an excited state, relaxes fast from 
this state, but the relaxation ends without depopulating 
completely the excited statei^. In these experiment the 
depopulation effects do not show any dependence on the 
resonance between the electronic inter-level energy sep- 
aration and optical phonon energy. These effects were 
recently interpreted to be a result of a non-geminate ex- 
citation of the electrons and holes in quantum dots^ii^ 
or as a manifestation of multiphonon scattering processes 
including acoustic phonons-^^. 

Besides the interpretations of the origin of these depop- 
ulation effects given in the experimental papersiiiSiiS, 
we propose the mechanism studied in this paper as an- 
other possible source of these effects. This mechanism 
assumes the absence of the phonon bottleneck effect: The 
formula Q for the relaxation rate does not contain the 
energy conservation delta function factor. It describes 
the electronic transitions between two electronic states. 



which are bound to a cloud of optical phonons. This co- 
herent phonons cloud is induced thanks to the electron- 
phonon coupling by the virtual transitions of the elec- 
tron between the two electronic states. Because these 
multiphonon states are not necessarily Fock states with 
well-defined number of LO phonons, the energy conserv- 
ing delta function does not appear in the relaxation rate 
formula. The coherent multiphonon states can be seen 
as classical macroscopic oscillations of the optical phonon 
modes. These macroscopic oscillations can be regarded 
as a source of an effective electronic Hamiltonian depend- 
ing explicitly on time, making the system nonconserva- 
tive. In such a system we do not speak about the energy 
conservation. This is the way how we can understand 
the electronic relaxation in the case of nonresonant sit- 
uation between the inter-level energy separation and the 
phonon energy, both in the case of the relaxation from 
the high energy level and in the case of the relaxation of 
the electron from the low-energy level to the high-energy 
one. 

The coupling of the electron to the phonon system via 
multiple scattering virtual transitions may cause a strong 
mixing of the two electronic states causes. Because of the 
effective lack of the energy conservation, the final state of 
the relaxation of the system need not to be identical with 
the low energy state, in the limit of zero temperature. 

Summing up, we have demonstrated the time depen- 
dence of electronic relaxation in the model quantum dot 
with two levels and a single electron. We demonstrate 
the incomplete depopulation of the excited state and we 
also demonstrate the appearence of the population of the 
excited state after the resonance excitation of the ground 
state. We present the characteristic times of these pro- 
cesses. These theoretical results are discussed in relation 
to the recent experiments and an additional interpreta- 
tion of these experiments is proposed. 
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